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, Abstract 

(N ■ 

^ , ' A formalism is introduced which may describe both standard linearized waves and grav- 

O [ itational waves in Isaacson's high-frequency limit. After emphasizing main differences be- 

■ tween the two approximation techniques we generalize the Isaacson method to non-vacuum 
spacetimes. Then we present three large explicit classes of solutions for high-frequency grav- 
itational waves in particular backgrounds. These involve non-expanding (plane, spherical 

I or hyperboloidal) , cylindrical, and expanding (spherical) waves propagating in various uni- 

T-H . verses which may contain a cosmological constant and electromagnetic field. Relations of 

^ I high-frequency gravitational perturbations of these types to corresponding exact radiative 

■ spacetimes are described. 

o : 
o ■ 
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O^' 1 Introduction 

u '. 

In classic work Isaacson presented a perturbation method which enables one to study prop- 

■ erties of high-frequency gravitational waves, together with their influence on the cosmological 
^ , background in which they propagate. It is this non- linear "back-reaction" effect on curvature of 

^ I the background spacetime which distinguishes the high-frequency approximation scheme from 

other perturbation methods such as the standard Einstein's linearization of gravitational field 
in flat space |2||^ or multipole expansions 01 that were developed to describe radiation from 
realistic astrophysical sources. 

The high-frequency perturbations were originally considered by Wheeler and then applied 
to investigation of gravitational geons by Brill and Hartle Isaacson's systematic study P 
stimulated further works in which his treatment was developed and also re-formulated in various 
formalisms. Choquet-Bruhat [7^IHI analyzed high-frequency gravitational radiation using a gen- 
eralized WKB "two-timing" method. Averaged Lagrangian technique which leads to Isaacson's 
results with less calculation was introduced by MacCallum and Taub [HIEI- Comparison of 
these approaches, and clarification of assumptions that have to be made in order to provide 
a consistent high-frequency approximation limit was also given by Araujo |lH ll2j. Elster |T3] 
proposed an alternative method that is based on expanding null-tetrad components of the Weyl 
tensor. Recently, Burnett developed a weak limit approach ^1] in which the high-frequency limit 
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can be introduced and studied in a mathematically rigorous way. These general methods have 
been, of course, applied to study explicit particular examples of high-frequency gravitational 
waves, see e.g. [HIHlinilEllini • 

On the other hand, many exact solutions of Einstein's equations are known which represent 
gravitational radiation. Among the most important classes are planar -waves |17(ll8j which 
belong to a large family of non-expanding radiative spacetimes |19l2Uj , cylindrical Einstein- Rosen 
waves [2^, expanding "spherical" waves of the Robinson- Trautman type |22|I23|. spacetimes 
with boost-rotation symmetry representing radiation generated by uniformly accelerated sources 
|24a25.26 , cosmological models of the Gowdy type j22j , and others — for comprehensive reviews 
containing also a number of references see, e.g., pH l lHn i lT H I29 | IH^ . 

However, there are only several works in which relation between exact gravitational waves 
and those obtained by perturbations of non-flat backgrounds has been explicitly investigated 
and clarified, see e.g. (10ti33ul6< . The purpose of our contribution is to help to fill this "gap". 

We first briefly summarize and generalize the Isaacson approach 1^ to admit non- vacuum 
backgrounds, the cosmological constant A in particular. Modification of Isaacson's formalism 
allows us to incorporate also standard linearized gravitational waves into the common formalism. 
Then, in section 3 we study properties of high-frequency gravitational waves in specific classes 
of spacetimes with special algebraic or geometric structure. In particular, we focus on waves 
which propagate in backgrounds with A 7^ 0. This is motivated not only theoretically but also 
by recent observations [HI] which seem to indicate that (effective) positive cosmological constant 
played a fundamental role in the early universe, but it is also important for its present and future 
dynamics. 



2 High-frequency approximation versus standard linearization 

Let us assume a formal decomposition of the spacetime metric g^i, into the background metric 
and its perturbation h^^, 

9fiu = lfj.u +£hfj_^ , (1) 

where, in a suitable coordinate system, 7^jy = 0(1) and h^u = 0(e) [by definition, / = 0{e'^) if 
there exists a constant C > such that |/| < Ce" as e 0]. The two distinct non- negative 
dimensionless parameters e and e have the following meaning: e is the usual amplitude parameter 
of weak gravitational perturbations whereas the frequency parameter e denotes the possible high- 
frequency character of radiation described by h^^. To be more specific, the parameter e <C 1 
characterizes (for e = 1) the amplitude of linearized gravitational waves in the ordinary weak 
field limit of Einstein's equations. The second independent parameter e = X/L represents, on the 
other hand, the ratio of a typical wavelength A of gravitational waves and the scale L on which 
the background curvature changes significantly. Isaacson's high-frequency approximation 
arises when A <C L, i.e. e <C 1 (and e = 1). Since L can be considered to have a finite value of 
order unity, we may write 0{e) = 0(A). 

To derive the dynamical field equations we start with the order-of-magnitude estimates 
which indicate how fast the metric components vary. Symbolically, the derivatives are of the 
order ~ l/L, dh ^ h/X, so that the following formulas 

7^, = 0(1) , V = 0{e) , 

lf.u,a = 0(1) , V,« = 0(1) , (2) 

l/iu,af3 = 0(1) , hf,^,a(S = 0(e"^) , 



2 



are valid. Next, we expand the Ricci tensor in powers of h, 



(3) 



where 
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The semicolons denote covariant differentiation with respect to the background metric j^i, , which 
is also used to raise or lower all indices. Considering relations @, the orders of the terms (@J) 
are 



Two limiting cases thus arise naturally. For the standard linearization (e <C 1, e = 1) the 
dominant term of Rfj_u{g) is R^p} = 0(1) which corresponds to the background 7^1, [to find, e.g., 
a vacuum spacetime metric g^y we solve R^p} (7) = 0] . Its first correction representing linearized 
(purely) gravitational waves is governed by 



which is a dynamical equation for perturbations h^y on the fixed background ^^y. The next term 
Rpv{l^ h) can then be used to define energy-momentum tensor of these gravitational waves, but 
the background metric is not assumed to be influenced by it. Improvements to this inconsistency 
can be obtained by iteration procedure. More rigorous but somewhat complicated solution to 
this problem was recently proposed by Efroimsky [HSj . 

In the high-frequency approximation (e <C 1, £ = 1) the dominant term is Rp'J = 0{e~^) 
which gives the wave equation © for the perturbations h^i, on the curved background 7^,^ 

(considering a vacuum full metric gpu)- The two terms of the order 0(1), namely Rp'u and 

(2) 

Rpu , are both used to give the Einstein equation for the background non-vacuum metric, which 
represents the essential influence of the high-frequency gravitational waves on the background. 
Of course, to obtain a consistent solution, one has to use both the wave equation and the Einstein 
equation for the background simultaneously. 

2.1 Linear approximation 

Interestingly, it follows that the wave equation for h^^, , which arises from the linear perturbation 
of the Ricci tensor in vacuum for both the above limiting cases e <C 1, e = 1, and e ^ 1, e = 1, 
is the same equation 0. In analogy with the well-known theory of massless spin-2 fields in flat 
space 12 we wish to impose two TT gauge conditions, 



= Oil), sR^^} = Oie^'e), e^^j = O(e^), e^R^ = 0{ee^). 



(5) 



R^'j{l,h) = 



(6) 



h 




(7) 
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In this gauge we arrive at the following wave equation 

V = - 24°i/3 ^'"'^ - <i ^''^ - ^% = , (9) 

where the operator <(> is the generalization of flat-space d'Alembertian. Contracting Q we 
obtain (/i'^^)'^.^ = 0, so that the condition (jSJ is always consistent with @. However, if we 
differentiate (^h^^, and use equations 0, ((2)), we find that 

{Oh,,r = - 2<;/3)^'^^ ' where (10) 

{Oh,.r = 0{e-^) , (R^X - 2<;/3)^'^^ = 0{e) . (11) 

Thus, in case of standard linearized waves (e = 1) there is an obvious inconsistency, except for 
backgrounds with a covariantly constant Ricci tensor (e.g., for the Einstein spaces). On the other 
hand, in the high-frequency limit (e = 1), the inconsistency between Q and is extremely 
small (the left and the right sides of 1)10(1 differ by where e <C 1). Moreover, for all background 
metrics of constant curvature the equations are fully consistent. This is an important advantage 
of the equation © containing also terms of non-dominant order (namely those proportional to 
the Riemann or Ricci tensors), if compared to other "simpler" wave equations (e.g., h^y'^ -p = 0) 
for which the left and right sides of ()1()|) generally differ by only two orders of magnitude. 



2.2 Nonlinear terms and the effective energy-momentum tensor 

Before considering the second-order terms we now extend the formalism to be applicable to a 
larger class of spacetimes with (possibly) non- vanishing energy-momentum tensor T^y. Namely, 
g^y need not be a vacuum metric (as only considered in P) but it satisfies Einstein's equations 

R^y{g)=M%y{g,ip) . (12) 

Here T^y = T^y — \g^j,yT^ such that T^y{g, ip) depends on non-gravitational fields ip and on the 
full metric gf^y but it does not contain the derivatives of g^y . Note that this admits as particular 
cases a presence of electromagnetic field, and also Einstein spaces when T^jy = -^Kg^y. Under 
the assumptions ((H) valid for the decomposition (P) we expand the equation (fT^ as 

<J(7) + eR^^Kl. h) + e24i)(7, h) + ...= (13) 
8^ [ (7, + e r(i) (7, h, ^) + fj^J (7, (^) + . . . ] , 

where T)il!^ (7, (/?) = f^y{'y,ip), and the remaining terms on the right-hand side are linear and 
quadratic in h, respectively. The orders of magnitude of the terms in the expansion of the Ricci 
tensor have been described above, cf. ©. For the energy- momentum tensor one obtains 

r(0) = O(l), f« = 0(6), f(l^ = 0{e'). (14) 

For ordinary linearization we thus get the equations R^^ = SvrT^"^ in each order n = 0, 1, 2, . . .. 
For the high-frequency approximation we obtain from (|13|) in the leading order 0(e~^) the 
equation © which is identical with the wave equation in the vacuum case. The second-order 
contributions, that are 0(1), represent an influence of the high-frequency gravitational waves 
and matter fields on the background, 

4°) (7) - SvtTW (7, ^) = (7, h) . (15) 
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This equation (which in case of a vacuum spacetime reduces to the Isaacson result) can be 
rewritten in the form of Einstein's equation for the background as 

G(0) (7) - Svr (7, ^) = - [41) (7, h) - i7M.«^'^ (7, h)] ^ 8vr T^J^ . (16) 

This defines the effective energy- momentum tensor T^J^ of high-frequency gravitational waves. 

2.3 Gravitational waves in the WKB approximation 

In the following we shall restrict ourselves to the Isaacson approximation (e = 1, e <C 1), i.e. on 
study of high-frequency gravitational waves on curved backgrounds. Inspired by the plane-wave 
solution in flat space, the form /i^^^ = ^e^,yexp(z0) of the solution is assumed. The amplitude 
^ = 0(e) is a slowly changing real function of position, the phase (/> is a real function with a large 
first derivative but no larger derivatives beyond, and e^i, is a normalized polarisation tensor field. 
The above assumption, introduced in pi, is called the WKB approximation, or the geometric 
optics limit |^ . The wave vector normal to surfaces of constant phase is = cf)^^ and the orders 

of various relevant quantities are R^^^^ = 0{1), A^i_i. = 0(e), = 0(e^^), and k^-,^ = 0(e~^). 
Substituting this into the conditions ((71), ©, and the wave equation Q we obtain, in the two 
highest orders which are gauge invariant, 

k^k^ = , k^e^i/ = , k e^jy-Q, = , 

e'^%, = l, 7A'-e^, = 0, (A^kP)^ = Q. (17) 

These express that a beam of high-frequency gravitational waves propagate along rays which are 
null geodesies with tangent /c'^, with parallelly transported polarization orthogonal to the rays. 
Moreover, using the WKB approximation of T^J^ and the Brill-Hartle averaging procedure 
(HI (which guarantees the gauge invariance) Isaacson obtained for gravitational waves in the 
geometric optics limit the energy-momentum tensor 

Tj^if = "^-^^n^u ■ (18) 

The energy- momentum tensor of high-frequency waves thus has the form of pure radiation. This 
fully agrees with results obtained by alternative techniques [5linill4j. 

3 Examples of high-frequency gravitational waves 

Now we present some explicit classes of high-frequency gravitational waves. These are obtained 
by the above described WKB approximation method considering specific families of background 
spacetimes with a privileged geometry. 

3.1 Non-expanding waves 

As the background we first consider the Kundt class |19|l2Sj of non-expanding, twist-free space- 
times in the form |36j 

ds^ = Fdv? -2-^dudv + —r (dx^ + dy2) , (19) 
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with 



P 

Q 



1 + 2 + y ) ' 



(20) 



(0% 

p2 



where a, /3, and e are constants (without loss of generahty e = or e = 1), Ci, C2 and D are 
arbitrary functions of the retarded time n, and H{x, y, u) is an arbitrary function of the spatial 
coordinates x, y, and of u. 

In particular, these are Petrov type N (or conformally flat) solutions of Einstein's equations 
with cosmological constant A when a = —j3 = gA and D = —2j3e + Cf + C2, see e.g. jSZlEEl 
I391l36j . Such metrics represent exact pure gravitational waves propagating along principal null 



direction if H satisfies the equation P {H^xx + H 



yy. 



+ IAH 



0. However, in our treatment 



here the function H does not describe exact gravitational waves but rather it characterizes the 
influence of high-frequency perturbations on the background metric, which is assumed to be 
initially given by (UH), ^ with H = 0. 

We consider the phase of high-frequency gravitational waves given by (/> = (f){u), and we seek 
solution in the WKB form, namely 



h 



exp 



(211 



where the amplitude A and polarization tensor e^^, are functions of the coordinates {u,v, x,y} . 
The corresponding wave vector is /c^ = ((/>, 0,0,0), where the dot denotes differentiation with 
respect to u. Applying now all the equations (fTTj) we obtain 
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(22) 



-'fj.u 



The fact that the amplitude A is independent of the coordinate v expresses non-expanding 
character of the waves. The special polarisation tensors, denoted as + and x, are analogous 
to those used in the standard theory of linearized waves in flat space. A general polarisation is 
easily obtained by considering e^,y = a + 6e^^, where a'^{u, x, y) + 6^(n, x, y) = 1. 

Using the Einstein tensor for the metric (|19j) with the cosmological term in equations (|16|) and 
()18() . we determine the reaction of the background on the presence of the above high-frequency 
gravitational perturbations, namely 



Q 
P 



dx'^ dy'^ 




H{u,x,y) = jA'^{u,x,y)(p^ 



(23) 



Notice that A = 0{e) and (f) = 0{e ^). Therefore, the influence of high-frequency gravitational 
waves on the background, represented by the function H, is of the order 0(1). These approximate 
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solutions can obviously be compared to specific exact radiative vacuum solutions which are given 
by H solving the field equation H23() with a vanishing right-hand side (when ^ = 0, i.e. high- 
frequency perturbation waves are absent). 

The above waves are non-expanding with the wave-fronts u = const, being two-dimensional 
spaces of constant curvature given by a = |A , cf. (|19|) . For A = these are plane- fronted waves, 
for A > they are spheres, and for A < hyperboloidal surfaces. 

Another interesting subclass of the Kundt spacetimes of the form (|19|) . (|20|) are explicit 
Petrov type // (or more special) metrics given by /? = a, e = 1, C = and D = 2(A — a), 
namely 



1 



2{A-a)v^-H du'' -2dudv + —r{dx^ + dy^) . (24) 
J 

For H = these are electrovacuum solutions with the geometry of a direct product of two 
2-spaces of constant curvature, in particular the Bertotti-Robinson, (anti-)Nariai or Plebahski- 
Hacyan spaces [TOU^EniEH] . see e.g. [SIISS]- Considering again (^1]) we obtain the results (|^ 
as in the previous case. However, the reaction of high-frequency waves on the background is now 
different. It is determined by the equations (|16() and ()18() with the energy- momentum tensor 
consisting of a cosmological term plus that of a uniform non-null electromagnetic field described 
by the complex self-dual Maxwell tensor F^'^ = 4<I>i(m[^?fi'^l — k^^^F^), where <I>i = — ^ e**^, 

c = const., and m = P d^, h = d^, 1 = dy + du form the null tetrad. Straightforward calcu- 
lation gives 

--lA\u,x,y)^^ . (25) 




This result is analogous to the equation (|23j) . but the present situation is now more compli- 
cated since the background spacetime is not vacuum but it contains electromagnetic field. (In 
fact, the term with the cosmological constant A in (|23j) has been entirely compensated by this.) 
Therefore, we have to analyze the perturbation of the complete Einstein-Maxwell system, and 
its consistency. 

The Einstein equations in the two highest orders © and (jlbj) have already been solved. 
We will now demonstrate that the Maxwell equations are also satisfied in the high-frequency 
limit, namely F^'^^,^ = 0(e), where | denotes the covariant derivative with respect to the full 
metric ^^j/. Indeed, using antisymmetry of F'^'^ we can write F^"\y = F^^ + ^g'^^ Qap^vF^" ■ 
Considering © and the gauge condition Q we obtain g°^^gap,v = 1°^^la(3,v — hf^^hap^u + O(e^) 
because j'^'^hap,u - h'^^jafs,^ = (^^13).^ - 2 h°''^'yai3;u = , so that 

F'^'i, = F'^^. - ^h'^^Kp^^Ff^" + 0(e2) . (26) 

Consequently, if the original background represents an electrovacuum spacetime, F'^'^-i, = 0, the 
Maxwell equations F^'^^^^ = 0{e) for the full metric are satisfied in the dominant order 0(1) in 
the high-frequency limit e ^ 1. In addition, the field equations are valid also in the next order 
0(e) for the new electromagnetic field 

^^- = (l+i/i°/3/i„^)F^- , (27) 

smce usmg 

(EHl) we obtain ^"""1^ = 0{e^). Starting fr om an electromagnetic field F^'^ satisfying 
F^^ -y = with respect to the background metric 7^1,, we can thus construct the electromagnetic 
field J-^'^ which satisfies the Maxwell equations J^^'^i^ = O(e^) with respect to the full metric g^y 
in the presence of high-frequency gravitational waves. Both the Einstein and Maxwell equations 
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are then satisfied in the two highest perturbative orders. Interestingly, these results hold for 
high-frequency perturbations of any "seed" electrovacuum background spacetimes. 

In particular, if the backgrounds are direct product spacetimes (|24|) for H = with uniform 
non-null electromagnetic field $i = const, then high-frequency gravitational waves (|21|) . ()22() 
introduce H which is given by equation (|25|) . According to (|27|) . the electromagnetic field is 
perturbed by the term proportional to h°'^haf3 = J^e^'^'^ = O(e^), see ((T7j) . namely 



$1 



1 + ^A^{u, X, y) e 



2i(j>(u) 



(28) 



This remains non-null but it is no longer uniform. The full spacetime thus describes non-uniform, 
non-null electromagnetic field plus the null field of high-frequency gravitational waves. 



3.2 Cylindrical waves 

Next we consider the class of cylindrical Einstein-Rosen waves. 



(29) 



If the functions '4'{t, p) and j{t, p) satisfy the corresponding field equations (see, e.g. [U, ^5] , 
or equations (jHSJ-llSSl) below) these are exact radiative spacetimes of the Petrov type I. We con- 
veniently define double null coordinates u = -^{t — p) and v = '^{t + p)] in these coordinates 
{u, V, if, z} the metric takes the form 



(30) 



We assume this to be the class of background universes into which we wish to introduce 



high-frequency gravitational waves. We assume again 



kj. 



(fT7|) we obtain 



,0,0,0), i.e. the WKB perturbation of the form ((21 



i){u) implying the wave vector 
By applying all the conditions 
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U{u) 

7^ 
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-2V' 
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(31) 



notice that v 



\ 

/ \ 

1 , , 

2^"-^) 1 

\ 1 / 

\/2/5 > 0. Thus the perturbative solution is given by 
U(u) 



h 



U 



exp 



(#(^)) . 



(32) 



The back-reaction on the background (contained in a specific modification of the metric 
functions 7 and -0) is given by the following equations, cf. ((TH|) . 



1{v 



(v - u) 7,„ = -^( 

{v - u) - = , 

{^,v - i',u) = . 



1^' 



1 



U 



(33) 
(34) 

(35) 
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Interestingly, this set of equations is consistent: by differentiating equation (|33jl with respect to 
V, equation (|5H) with respect to u, and combining them, one obtains H35() provided the ampUtude 
A{u, v) satisfies the equation 

-u)A^) =0. (36) 



However, this is automaticaUy satisfied for the amphtude (|31j) . It is thus quite simple to intro- 
duce gravitational waves in the WKB approximation into the cylindrical spacetimes pOf) . If the 
functions 7 and ip representing the background are solutions of the vacuum equations [i.e. (|38|l - 
()35() with a vanishing right-hand side of H33() ] then for introducing high-frequency gravitational 
waves it is sufficient just to alter the function 7 as 

7(u, v) 7(u, v) + j{u) , (37) 

where 



du 16 



In particular, when ■0 = = 7 the background ()29|) is a flat Minkowski space. By assuming non- 
trivial 7 we obtain Petrov type spacetime with high-frequency gravitational waves which have 
cylindrical wave- fronts. In a general case this perturbation is propagating in the background 
which is the Einstein- Rosen cylindrical wave of Petrov type I. The effect on background is given 
by the relation (|n5|) where U{u) = 0(e) is an arbitrary amplitude function. 

The above described perturbations depend on the null "retarded" coordinate u so that the 
high-frequency gravitational waves are outgoing (p is growing with t, on a fixed u). However, 
since the background metric H3U|) is invariant with respect to interchanging u with v, it is 
straightforward to consider also ingoing perturbations by assuming the phase to depend on the 
"advanced coordinate" v, namely 

hi,u = e^iv exp («(/>(«)) . (39) 

^Ju — V ^ ' 

Then the term proportional to A?(fP' will appear on the right-hand side of equation (|34)) instead 
of This results in an interesting possibility to introduce ingoing high-frequency gravitational 
cylindrical waves into the background of outgoing Einstein-Rosen waves just by assuming 7(11) 
in (|S7|) such that 

^^'''^ +i^V202^ (40) 



dv 16 



or vice versa. 

Moreover, all the above results can further be extended to a class of generalized Einstein- 
Rosen (diagonal) metrics [^115] which describe G2 inhomogeneous cosmological models, 

ds' = e2^-2V'(-dt2 + d/>2) + e^^dz^ + t\-^^d^^ . (41) 

If the three-dimensional spacelike hypersurfaces are compact, the corresponding model is the 
famous Gowdy universe with the topology of three-torus j27|l29) . In the double null coordinates 
just one component of the metric is now different from H30() . namely g^p^p = ^(v-\-u)'^e~'^'^^'^''"\ The 
only modification of the above results (in the double null coordinates) consists of replacing the 
factor {v — u) with {v-\-u), and each derivative with respect to u changing sign (e.g. 7^^ —l,u 
or "ip^uv —'ip,uv)- High-frequency gravitational waves in inhomogeneous cosmologies of the form 
(|^T|) can thus easily be constructed. 
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3.3 Expanding waves 



Finally, we assume that the background is an expanding Robinson- Trautman spacetime. The 
metric (generally of the Petrov type II) in the standard coordinates has the form, see e.g. 

isaEsiEHiisn], 



K -2r(lnV)^u-2 



m 



A 



2dudr 



p2 



(d7?2 + d^^ 



(42) 



where K = A(lnP), A 



+ 



, and m{u). When V{u,r],^) satisfies the Robinson- 
4:^,u = 0, the metric (|42|) is an exact vacuum solution 



Trautman equation AK + 12 m (In V)^u 
of the Einstein equations. 

In view of the existence of privileged congruence of null geodesies generated by dr we intro- 
duce the phase cp = 4>{'^) and the wave vector = (0, 0, 0, 0) of high-frequency gravitational 
waves. We again assume the WKB form (|2T|) of the solution. Applying the equations ([T7jl we 
obtain 



A 
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71^ 
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\ 1 y 



(43) 



A general solution has the form h^y = r -"^[/(tt, ry, ^) e^,/ exp ^i(/)(u)^, where U{u,r],^) and (j){u) 
are arbitrary functions, and Cf^u = ae^^ + be^y with a^{u,ri,^) + }p'{u,ri,£) = 1. Introducing the 



amplitudes = aU, = bU for both polarizations, we can write the solution as 



h 



exp 



(#(n)) . 



(44) 



If the wave-surfaces r = const. ,u = const, with the metric dl"^ = V^'^{d'rf' +d^^) are homeomorfic 
to 5^, the waves can be interpreted as "spherical". In the asymptotic region r ^ oo such 
solutions locally approach plane waves |16j . 

The reaction of the waves on background is determined by the equations (jlOf) and H18|) with 



(0) 



-^A7^i/. From the only nontrivial component we immediately obtain the following 



equation 



^ + 3m (lnP),„ + 1 = ^ [(C/+)2 + ([/ 



x\2 



(45) 



where m{u), (p{u), whereas the remaining functions depend on coordinates {u, tj, (^}. Notice that 
this is independent of the cosmological constant A. 

The expressions (|44|) , (|45p agree with results obtained by MacCallum and Taub 01 or recently 
by Hogan and Futamasc 16 who used Burnett's technique [^. Our results, which were derived 
by a straightforward approach, are slightly more general because they are not restricted to a 
constant frequency cj) = const. Particular subcase of the Vaidya metric has already been studied 
before by Isaacson 1 and elsewhere [H]. 
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4 Conclusions 



The Isaacson approach to study high-frequency perturbations of Einstein's equations was briefly 
reviewed and compared with the standard weak- field limit. In our contribution we generalized 
Isaacson's method to include non-vacuum spacetimes, in particular an electromagnetic field 
and/or a non- vanishing value of the cosmological constant A. Then we explicitly analyzed pos- 
sible high-frequency gravitational waves in three large families of background universes, namely 
non-expanding spacetimes of the Kundt type, cylindrical Einstein-Rosen waves and related in- 
homogeneous cosmological models (such as the Gowdy universe), and the Robinson- Trautman 
expanding spacetimes. These backgrounds are of various Petrov types. For example, high- 
frequency gravitational waves can be introduced into electrovacuum conformally flat Bertotti- 
Robinson space, type D Nariai and Plebahski-Hacyan spaces, their type and type // gener- 
alizations, or into algebraically general Einstein-Rosen universes. 

For construction of high-frequency gravitational perturbations we have employed the fact 
that all these spacetimes admit a non-twisting congruence of null geodesies. The corresponding 
tangent vectors k^^ are hypersurface orthogonal so that there exists a phase function </> which 
satisfies (j)^^ = k^. The last equation in (|T7j) can be put into the form Jj(ln^) = —0, where I is 
the affine parameter, and = \k^^^ is the expansion of the null congruence. This determines the 
behaviour of the amplitude A in the above spacetimes ()22() . H43() . The remaining equations 
()17() enables one to deduce the polarization tensors. 

It has been also crucial that all the classes of spacetimes discussed admit exact solutions 
with the energy-momentum tensor of pure radiation, i.e., G^u — Svr T^u = | J?k^ky, where T^y is 
either constant (representing the cosmological constant) or it describes an electromagnetic field. 
The relation between high-frequency perturbations and exact radiative solutions of Einstein's 
equations in each class is thus natural. In particular, it is possible to determine explicitly the 
reaction of the background on the presence of high-frequency gravitational waves. 
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